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ABSTRACT 

The globular cluster sy stem of a typical spheroidal galaxy makes up about 0.25% of the total galaxy 
mass (|McLaughhi]l ll999r This is roughly the same mass fraction as contained in the nuclear star clus- 
ter (or stellar nucleus) present in most nearby low-mass galaxies. Motivated by this "coincidence" , 
this Letter discusses a scenario in which globular clusters of present-day galaxies are the surviving 
nuclei of the dwarf galaxies that - according to the hierarchical merging paradigm of galaxy forma- 
tion - constitute the "building bl ocks" of present-day massive galaxies. This scenario, which was 
first suggested bv iFreemanI |1993'), has become more attractive recently in the light of studies that 
demonstrate a complex star formation history in a number of massive globular clusters. 
Subject headings: globular clusters: general 



1. INTRODUCTION 

The origin of globular cluster (GC) systems of present- 
day galaxies continues to provide a rich topic for investi- 
gation, speculation, and debate. Enabled mostly by the 
spatial resolution and precision photometry of the Hub- 
ble Space Telescope (HST), recent studies have estab- 
lished that GCs are not uniform in their properties, but 
display a rather wide range in color and metallicity. The 
color distribution of GCs in early-type galaxies nearly al- 
ways shows a doublc-peakcd structure ( Peng ct al. 2006) , 
motivating a commonly accepted view of GC systems 
as "bimodal" . The situation in disk-dominated galax- 
ies is less clear, partly because the analysis is limited 
by the smaller number of GCs in less massive galaxies. 
While both the Milky Way and M31 show clear signs 
of a double-peaked color distribution, the GC popula- 
tions in more distant spirals are often divided into "blue" 
and "red" subpopul ations hy definition through color se- 
lectio n criteria (e.g. iGoudfrooii et al.ll2003l : iRhode et al.l 
l2007f ). Whether or not "bimodality" is universal, the 
broad range in GC color makes it clear that GCs within 
a given galaxy have formed neither coevally nor under 
the same physical conditions. 

Moreover, the ongoing formation of young massive 
cluster ( YMCs) observed in nearby starbursts such 
as M82 (iMcCradv fc GrlOmJ [2007h or NGC 4038/9 
(|Mengel et al.l |2002[ ). as well as the presence of some 
intermediate-age GCs in, e.g., t he Magellanic Cloud s 
(IGoudfrooii et al.l[200fih or M33 (jChandar et al.l[200l . 
suggest that GC formation was not confined to the early 
universe, and possibly continues until the present day. 
Although it seems plausible that individual YMCs can 
evolve into GC-like structures, it is less clear that entire 
GC systems can be created in this way. It may well be 
that nature has more than one metho d to make GCs. 

The recent, comprehensive review of iBrodie fc Stradeii 
(|200Q) summarizes the proposed scenarios for GC for- 
mation. The most widely used are i) in s itu formation 
throu gh multi-phase dissipational collapse ([Forbes et al.l 
Il997f ). ii) massive star formation trigge red by major 
disk-disk mergers (lAshman fc Zepj|1992[ ). and iii) dis- 
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sipationless accretion of GC systems during the assem- 
bly of massive ga laxies through mergers at high redshift 
(jCote et al.lll998D . 

These models all have considerable problems in fully 
explaining the observed properties of GCs. For exam- 
ple, the multi-phase dissipational collapse model invokes 
two distinct phases of GC formation to explain color bi- 
modality, but cannot readily provide a physical mech- 
anism that would explain the dormant period between 
them. The major merger model, on the other hand, fails 
to account for the lack of intermediate-age GC systems 
and the systematically higher specific frequency in gE 
galaxies which are presumed to be the end products of 
major disk-disk mergers. Fi nally, the dissipationless ac- 
cretion scenario (|Cote et al.i ri998). as well as more recent 
studies emplo ying semi-analytic models of galax y forma- 
tion (e.g. Beas lev et aLll2002t iRhode et al.ll2005D . appear 
able to explain the observed color distributions of GC 
systems. However, these models do not explain how the 
GCs of the first galaxies formed, and thus merely shift 
the problem to earlier times. 

Most importantly, none of these scenarios can explain 
the observed complexity in the star formation history 
of some GCs. While most "typical" GCs studied so far 
show no obvious signs of multiple stellar populations (e.g. 
ISarajedini et al.ll2007f l. the situation is different for the 
most massive GCs. For example, the high mass and mul- 
tiple stellar populations of lo Cen have long triggered 
speculation that it is the remnant nucleus of a satel- 
lite galaxy which was destroyed when merging with the 
Milky Way (jNorris et all 119971 : iLee et al.|[l999l ). More- 
over, recent observations have established that oj Cen is 
not a unique case, and that the high-mass end of the 
Milky Way GC system commonly s hows evidence for 
multiple stellar populati ons (Moehle r fc Sweigartj 120061 : 
iPiotto et al.l l2007t ICaloi fc D'Antonall2007t ). These re- 
sults have challenged the traditional view that all GCs 
are simple stellar populations born in a single, short-lived 
( "instantaneous" ) burst of star formation. 

This Letter revisits and expands on an alternative for- 
mation s cenario for GCs which was first suggested by 
iFreemaiil ((1993) . The underlying idea is that GCs form as 
the nuclei of dwarf galaxies in the early universe, and are 
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accreted over time as their hosts merge into larger struc- 
tures. The paper is s tructu red as follows: § [2] briefly de- 
scribes the iFreemarJ (|1993D scenario, and mentions some 
recent observational developments that are relevant in 
this context. § [3] reviews what is known about the for- 
mation of NCs in late-type disk galaxies. § [4] then dis- 
cusses some quantitative aspects of the proposed scenario 
and presents a first order plausibility check. Finally, § [5] 
concludes with some open issues and predictions. 

2. GLOBULAR CLUSTERS AS DWARF NUCLEI 

Based on the observed lack of an abundance gradi- 
ent in the GC syste m of the Milky Way outer halo, 
ISearle fc Zi^ ^ME) suggested that GCs originate in 
small "fragments" (i.e. proto-galaxies) which are ac- 
creted over an extended period of time into a larger host. 
However, they did not discuss the specifics of where and 
how the GC-precursors form inside the fragments. 

Building on this scenario, iFreemai] ((1993) suggested 
that nucleated dwarf galaxies and their compact stellar 
cores might be the present-day analogs of those Searle- 
Zinn fragments that escaped disruption in a merger 
event. He also pointed out that chemical inhomo- 
geneities in present-day GCs might be explained by self- 
enrichment due to the strong stellar winds observed in 
starbursting galaxy nuclei. The general idea behind this 
scenario has recently received increased attention be- 
cause it might provide a natural explanation for the ob- 
served multiple stellar populations in some GCs, espe- 
cially when taking into account two other recent obser- 
vational results: 

1) stellar nucleation is not limited to dE galaxies, but 
is also common in low-m ass disk galaxies, as shown by re- 
cent surveys with HST (jCarollo et al.lll997l : iBoker et all 
|2002| ). These nuclear star clusters (NCs) are structurally 
very similar to GCs, and nearly always contain multiple 
stellar populations. They are therefore attractive candi- 
dates for the precursors of at least those GCs that have 
complex star formation histories. In fact, it is difficult 
to explain how a stellar cluster can accrete (or retain) 
enough gas for multiple star formation episodes if it is 
not located at the bottom of a potential well such as the 
nucleus of a galaxy. 

2) the mass of a GC syste m appears to be pr oportional 
to that of its host galaxy (jMcLaughlinI Il99"9h . As dis- 
cussed further in § 21 this is expected if GCs indeed form 
as the stellar nuclei of "building blocks" which merge into 
massi ve present - day g alaxies. 

The iFreemaJ (|1993h scenario can thus be modified as 
follows: present-day GCs have formed as NCs in the cen- 
ters of low-mass disk galaxies in the early universe where 
they grew in mass through a series of discrete star forma- 
tion events as long as the disk structure of their host re- 
mained undisturbed by any merger activity. This growth 
process ends when the host galaxy disk is disrupted in a 
merger or tidal interaction. Due to its large stellar den- 
sity, the (former) NC survives the merger, but because it 
is no longer located in a potential well, it experiences no 
further gas inflow. No further star formation takes place 
in the NC, and it will simply age passively as a GC in the 
merger product. It should be noted explicitly that the 
presence of multiple stellar populations is not a necessary 
consequence of this process: if a NC was accreted early 
enough, it might have experienced only a single star for- 



mation episode and thus contain only a single generation 
of stars. 

3. NUCLEAR STAR CLUSTERS 

Over the last decade, a number of studies - both 
via imaging and spectroscopic observations - have con- 
tributed to the following picture of present-day NCs: 

1) NCs arc common: the fraction of galaxies with an 
unambi guous NC detectio n is 75% in late-type (Scd-Sm) 



spirals llBoker et al.ll2002D. 58% in earlier- type (SO-Sbc) 
spirals (jBalcells. Graham, fc Peletieill2007D. an d 70% in 
spheroidal (E & SO) galaxies (lCote_eLaUl2006D^ 

2) The structural properties of NCs in late-type disks 
are indistinguishable from those of GCs, and thus cer- 
tainly do not rule out an evoluti onary connectiori : typi- 
cal half-light radii are 2 - 5pc ( Boker et al.ll2004f) . The 
size distribution of NCs in spheroidal galaxies also shows 
a large overlap with that of Milky Way GCs, although 
some present- day (d)E nuclei are clearly larger th an the 
"typical" GC (|Geha et al.ll200l iCote et al.ll2006n . 

3) Despite their compactness, present-day NCs are 
rather n iassive: their typica l dynamical mass is 10^ — 
10^ Mq (jWalcher et al.ll2005D which is at the high end of 
the GC mass function. Because of their very high stel- 
lar densities, there is little doubt that NCs are difficult 
to disrupt, and that they will remain structurally intact 
even after the disruption of their host galaxy. 

4) The star formation history of NCs in disk galaxies is 
complex, as evidenced by the fact that they have stellar 
populations compris e d of multiple gen erations of stars 
(jWalcher et an[200fil : iRossa et aLl[20M) . In the NCs of 
late-type disks, the youngest stellar generation is nearly 
always younger than lOOMyr which explains why present- 
day NCs are much more luminous than GCs (Boker et aL] 
120021 : ICote et all (20061 ). Nevertheless, NCs always con- 
tain an old (~ Gyrs) stellar population, suggesting that 
they have been in place for a long time. 

Taken together, these facts suggest that NCs have built 
up the bulk of their presen t-day mass through a series 
of star formation episodes ()Walcher et al.l |2006[ ) . This 
repetitive "rejuvenation" can only occur as long as gas is 
funneled towards the nucleus, and thus no longer takes 
place after the NC has been removed from the galaxy 
nucleus or after the gas supply has otherwise been dis- 
rupted, e.g. through transformation of a late-type disk 
into a dE galaxy. 

Three recent and independent studies of NCs in dif- 
fereii t galaxy ty pes fRo ssa et al.ll2006l : IWehner fc Harris! 
120061 : ICote et a l. 2006) have established that NCs obey 
similar scaling relationships with host galaxy proper- 
ties as do supermas sive black holes. This has led 
iFerrarese et al.l ()2006l ) to suggest that NCs and black 
holes might constitute two different incarnations of a 
"compact massive object" , the formation of which is in- 
herent to galaxy assembly, and which contai ns ^ 0.2% 
of the host galaxy (bulge) mass. Although Ba lcells et al.l 
(|2007t l have questioned the simple linearity, they confirm 
that any NC typically contains 0.1% — 0.3% of the bulge 

^ All these numbers are likely lower limits, albeit for different 
reasons. In the latest-type disks, it is sometimes not trivial to 
locate the galaxy center unambiguously so that no particular source 
can be identified with it. In contrast, many early-type galaxies have 
very steep surface brightness profiles (SBPs) that make it difficult 
to detect even luminous clusters against this bright background. 
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mass. 

A possible theoretical explanation for a similar pro- 
portionality also in trul y bulgeless disk galaxies has re - 
cently been offered by lEmsellem fc van de VenI (|2007h 
who point out the compressive nature of tidal forces in 
a shallow gravitational potential, e.g. in an exponential 
disk. Based on their simulations, a NC formed through 
this mechanism should contain ~ 0.1 — 0.5% of the host 
galaxy mass. In this model, matter infall onto the nu- 
cleus occurs naturally in bulge-less disks, possibly ex- 
plaining why NCs with recent star formation have only 
been found in in late-type spirals. The following section 
investigates the possibility that many, if not all, GCs are 
the surviving nuclei of dwarf galaxies that were "swal- 
lowed" during the assembly of present-day galaxies. 

4. A SIMPLISTIC SCENARIO 

A prerequisite for the proposed scenario is that NCs 
are present already in the first galaxies, i.e. the "build- 
ing blocks" for galaxy assembl y. This is ii ot an implausi- 
ble proposition. For example. ICenI (|2001h has suggested 
that the reionization of the universe produces inward, 
convergent shocks in dark matter "mini-halos" . These 
shocks compress the baryonic gas which in turn becomes 
self-gravitating and undergoes star formation, forming a 
dense stellar cluster in the process. These "seed clus- 
ters" naturally define a ceii t er of gravity for the assem- 
bly of a gas disk. As ICenI ()200lD points out, the NCs 
of high-redshift mini-halos formed through this process 
are plausib le progenitors of present-day GCs (see also 
iMoore et al . 2006). 

The subsequent evolution of these "primordial" NCs 
in their parent halos an d the expected tidal fields has 
been studied in detail bv iMashchenko fc Sill^ ()2005al lbl). 
These authors made the a priori assumption that 0.88% 
of the total halo mass is concentrated in a NC, a some- 
what higher fraction than what has been found for 
prese nt-day NCs . These authors conclude, in agreement 
with ICenI ()200lD , that the predicted properties of such 
NCs match those of present-day GCs, independent of 
the exact value of the initial NC mass fraction. 

They also point out that current observations are in- 
sufficient to rule out significant amounts of dark matter 
in GCs which might be expected if GCs indeed form at 
the centers of mini-halos. So far, however, dark mat- 
ter has not been convincingly detected in GCs, although 
th e fiat velocity dispe rsion profile of some GCs found 
bv lScarpa et al.l (|2007l ) appears consistent with the pres- 
ence of dark matter (or, as those authors prefer, with 
modified Newtonian dynamics). In any case, even the 
apparent absence of dark matt er in GCs does n 't fully 
invalidate the idea, because IMashchenko fc Silli (j2005bl ) 
show that, depending on the density profile of the halo, 
dark matter can be lost almost completely by tidal strip- 
ping. 

In this paper, thelci^ ((200lh scenario is modified and 
extended by allowing for the expected secular evolution 
of NCs that occurs as long as they remain in the center 
of their parent halo. This secular evolution offers a nat- 
ural explanation for the diversity in stellar populations 
within and across GCs. For a first-order assessement, the 
following simplifying assumptions are made: 

a) all building block galaxies are identical in that 
their baryonic matter is organized in a gas-rich disk of 



~ 10^ — 10^ Mq. This mass range is plausible because 
it makes dynamical friction an efficient mechanism to 
m erge small galaxies into a larger halo (see Eq. 7-27 
in iBinnev fc Tremaind Il987f) . It is also consistent wit h 
models of structure formation (e.g. iMoore et al]l2006[ ). 
as well as with obs erved masses of d warf galaxies in the 
local universe fe.g. lGeha et~all[200l . 

b) the center of each progenitor disk is occupied by a 
NC with ~0.2% of the host mass. As long as the disk 
escapes merger activity, its NC grows in mass via a series 
of nuclear starbursts fueled by gas flow from the disk onto 
the nucleus, similar to those observed i n nearby galaxies 
(e.g. iSchinnerer eran[2003l . [M)fiL[200l . 

c) once the disk experiences a destructive merger, the 
NC is captured in the halo of the merger product. From 
then on, it ages passively as a GC. Therefore, the later 
a galaxy is accreted, the more massive the newly added 
GC, and the more complex its star formation history. 

d) stars and gas of the merging disk will be incorpo- 
rated into the merger products' halo and disk, respec- 
tively, so that the total mass of the merger product is 
simply the sum of all the constituents. 

e) the potential for cluster formation induced by the 
merging process is ignored for the moment, i.e. GCs are 
only added through accretion of NCs. 

In this simplistic scenario, the total number of GCs 
is equal to the number of accreted building blocks. Be- 
cause any intermediate merger products are created from 
identical building blocks as described in assumption a), 
the number of GCs simply scales linearly with assem- 
bled galaxy mass, regardless of the details of the merger 
tree. This offers a natural explanation for the empir- 
ical result that all but the latest-type spirals (see SjS]) 
show a c onstant specific freque ncy of GCs: the observed 
T-value (Zepf & Ashmau'igga") is - 2 GCs per 10^ Mq 
([Goudfrooii et al. 2003; Rhode et al. 2007). This value is 
likely a lower limit because i) not all GCs will survive and 
ii) the faintest GCs might be missed in surveys. Taken 
at face value, however, it follows that the building blocks 
should have masses of < 5 x 10^ Mq, entirely consistent 
with assumption a) above. 

The primordial NC of such a galaxy should then con- 
tain 0.002 X 5 • 10^ Mq, i.e. < 10^ Mq, in good agree- 
ment with the smalle st observed NC mass of 8 x 10^ Mq 
(jWalcher et al.ll2005[ ). This apparent minim um NC mass 
is unl ikely to be a selection effect, because iBoker et al.l 
(|2002f ) found a faint-end cutoff in the NC luminosity 
function in a much larger, unbiased sample of NCs. Note 
that NC masses in this range lie well above the critical 
mass ~ 2 x 10^ Mq at which the GC mass function 
changes slo pe because less mas sive clusters are prone to 
disruption (|Fall fc Zhandl2b01h . It is therefore reason- 
able to expect that any such clusters will have survived 
until the present day. 

Nevertheless, depending on the details of their orbit 
and initial density profile, the captured clusters will expe- 
rience two-body relaxation and/or tidal stripping. These 
effects can easily red uce the cluster mass by a f actor of 2- 
5 over lOGyr (e.g. lVesperini fc Heggidl 19971) . An even 
highe r mass loss has been reported in some Milky Wa y 
GCs (Ide Marchi et al.l[2006l : Ide Marchi fc Pulonell2007D . 
making dynamical evolution a likely explanation for the 
fact that many GCs today have m asses below 10^ Mq 
(|McLaughlin fc van der MareH[2005f) . 
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If present-day GC systems are indeed made from NCs 
accumulated throughout the merging history of a galaxy, 
then the mass of any GC system should - to first or- 
der - contain the same mass fraction as each individual 
NO (because mass is conserved according to assumption 
c) above), i.e. ~0.2% of the total host galaxy mass. 
Intriguingly , this appea r s enti rely consistent with ob- 
servations: iMcLaughlinI (|1999D concludes that the for- 
mation efficiency of GCs is universal, at least for el- 
lipticals and the spheroidal component of disk galax- 
ies (inch the halo of the Milky Way), with a value of 
eel = MGCs/(Mgas + Af.tars) ~ 0.0026. It should be 
noted, however, that the situation is less clear when en- 
tire disk galaxi es (i.e. including the g as-rich thin disk) 
are considered (|Brodie fc StradeJl2006|) . 

Assuming for the moment that ed is indeed universal, 
and had the same value at all redshifts, the argument can 
be turned around: the putative GC system of a 10^ Mq 
building block galaxy should then contain a mass of only 
5 X 10^ M0. Unless this mass is concentrated in very few 
(one or two) clusters, none of the GCs in that system will 
have survived until the present day. Note that this also 
represents a challenge for the dissipationless accretion 
scenario, because in order to produce significant numbers 
of stellar clusters in the halos of early-universe galaxies 
with enough mass to survive until z—0, the value of ed 
must have been significantly higher in the past. 

5. OPEN ISSUES AND PREDICTIONS 

The toy model described in the last section is clearly 
simplistic, and there are a few points that require fur- 
ther investigation. For example, it has been assumed 
throughout the paper that the mass fraction of NCs re- 
mains constant over time at 0.2% of the host galaxy, 
as observed in present-day galaxies. For this assump- 
tion to hold, the gas accretion rate onto the host galaxy 
must match the mass growth of the NC. To first order, 
such a link does not appear implausible because there 
are a number of mechanisms that transport gas from the 
outskirts of a galaxy towards the center. However, it is 
presently unclear whether galaxy evolution models can 
indeed explain such "fine tuning" . 

Moreover, there is no doubt that the stellar popula- 
tions of GCs are generally old. For example, the Milky 
Way GCs have a ges of > 8 Gyr with a relatively small 
spread of < 3 Gyr ( De Angeli et al.ll2005[ ). Because in the 
proposed scenario, the age of the youngest stellar popula- 
tion in a GC denotes the time when its former host was 
disrupted, this implies that the merging of disk galaxies 
into the Milky Way was mostly complete 8 Gyr ago. Any 
NCs accreted after that time must have been residing in 
dE-type galaxies, where they had already been aging pas- 
sively for some time because there is no gas supplied to 
the nucleus in these galaxies. While this conclusion does 
not seem to be contradicted by any current observations, 
more detailed simulations of the Milky Way formation 
are required to test its plausibility. 

A third issue that is not explained by the proposed 
scenario is the double-peaked color distribution of many 
GC systems. In massive galaxies, where this effect is 
most pronounced, it may be possible to invoke additional 
cluster formation during the merger of two massive spi- 
rals (each of which has formed its own GC system ac- 
cording to the proposed scenario), but not in any mi- 



nor mergers that led to the build-up of the disk galaxies 
themselves. That star formation is indeed enhanced in 
major mergers is undisputed an d has been observed in 
a number of nearby objects (e.g. IWhitmore et al.l [19931 : 
IWhitmore fc Schweized ll995V However, the fact that 
there are hardly any GCs known with intermediate ages 
(few Gyr), requires that the majority of these mergers 
were completed by z ~ 2, possibly a challenge to mod- 
els of galaxy formation. On the other hand, it may not 
even be necessary to invoke the presence of two GC sub- 
populations with intrin sically different metallic ities to ex- 
plain color bimodality (|Yoon. Yi. fc Leell2006f ). 

One prediction that the proposed scenario can be 
tested against is that "pure" disk galaxies should have 
a systematically lower T- value. This is because they 
apparently had an uneventful merger history which en- 
abled them to avoid the formation of a bulge. Be- 
cause these galaxies likely assembled a larger fraction 
of their mass through slow accretion of intergalactic gas 
rather than via accretion of satellites, they should have 
gathered proportionally fewer NCs. Generally speak- 
ing, little is known about the GC systems of the latest- 
type spirals, although a few studies have been under- 
take n which seem to confirm this p rediction. For exam- 
ple, Ivan denBergh & Harris^ ()1982[ ) report the absence 
of GCs in the edge-on spiral NGC 891 which, however, 
might not be a good example of a bulge-less system since 
it is classified as Sb. Moreover. lOlsen et al.l (|2004f ) find 
only a handful of confirmed GCs in most late-type spi- 
rals of the Sculptor group, and report that the specific 
frequency of GCs is lower in late-type spirals than in 
bulge-dominated galaxies. Unfortunately, these authors 
do not present T- values which is a more robust diagnos- 
tic because of potential color variations between galaxy 
types. 

It is still uncertain when exactly NCs formed in the uni- 
verse and what the duty cycle of rejuvenation in those 
early nuclei was. In principle, one would expect that 
a large fraction of GCs show some level of complexity 
in their stellar populations. However, the earlier a NC 
was acquired, the fewer stellar generations it should con- 
tain, and the smaller the expected age differences be- 
tween those populations are. As mentioned in §[2l single- 
population GCs are certainly expected, especially if there 
was a co nsiderable delay between the original collapse of 
the ICenI {2001) model, and the onset of "rejuvenation" 
of the NC. 

On average, however, more massive GCs should have 
been accreted more recently than low-mass ones. Be- 
cause they spent a longer time in a galaxy nucleus, they 
experienced a higher number of star formation events, 
and contain both younger and more complex stellar pop- 
ulations. The fact that multiple stellar populations so far 
have mostly been detected in rather massive GCs such 
as NGC 2808 and uj Cen thus is not surprising. 

In summary, the speculative scenario for GC forma- 
tion laid out in this letter does not seem implausible, 
and offers a natural explanation for the observed mass of 
GC systems (Mqcs ~ 0.0025 • Afgai) and the presence of 
multiple stellar populations in massive GCs. 
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